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Abstract 

The semi-inclusive reaction e + p — > e + Xp was studied with the ZEUS detector at 
HERA using an integrated luminosity of 12.8 pb _1 . The final-state proton, which 
was detected with the ZEUS leading proton spectrometer, carried a large fraction 
of the incoming proton energy, xl > 0.32, and its transverse momentum squared 
satisfied pf, < 0.5 GeV 2 ; the exchanged photon virtuality, Q 2 , was greater than 3 
GeV 2 and the range of the masses of the photon-proton system was 45 < W < 225 
GeV. The leading proton production cross section and rates are presented as a 
function of xl, Q 2 and the Bjorken scaling variable, x. 
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1 Introduction 



Hadron-hadron collisions predominantly give rise to leading particles of the same type as 
those in the incoming beams and carrying a large fraction of the momentum of incoming 
particles. The spectrum of leading particles approximately scales with the centre-of- 
mass energy, a property known as limiting fragmentation [1]. The properties of the 
accompanying final state are also universal when studied as a function of the centre-of- 
mass energy available after excluding the leading particles [2,3]. 

Events with a final-state proton carrying a large fraction of the available energy, Xl, but 
a small transverse momentum, p?, have been studied in detail in high-energy hadron- 
proton collisions [4-6]. More recently, the HERA experiments reported measurements 
of the production of leading protons in ep collisions [7, 8]. Several mechanisms have 
been suggested to explain the production of leading protons. None of them are, as yet, 
amenable to calculations based on perturbative quantum chromodynamics (pQCD). This 
is, in part, a consequence of the small values of px of the leading proton which necessi- 
tates a non-perturbative approach. Some models [9-13] are based on the Regge formalism, 
with leading proton production occurring through t-channel exchanges, both isoscalar and 
isovector, notably of the Pomeron, pion and Reggeon trajectories. These exchanges medi- 
ate the interaction between the proton and the hadronic fluctuations of the virtual photon. 
Other models retain quarks and gluons as fundamental entities, but add non-perturbative 
elements, such as soft-colour interactions (SCI) [14]. Alternatively, the concept of fracture 
functions [15] offers a QCD framework in which to describe the leading baryon momentum 
spectra. 

This paper presents measurements of leading proton production in e + p collisions, e + p — > 
e + Xp, with a four- fold increase in statistics compared to an earlier measurement [7]. High- 
energy protons with low transverse momentum carrying at least a fraction xl=0.32 of the 
incoming-proton momentum were measured in the ZEUS leading proton spectrometer 
(LPS) [16]. The cross sections are presented as a function of the proton variables xl 
and p\. The dependence on the Bjorken variable, x, and on the photon virtuality, Q 2 , 
was also studied and compared to that of the inclusive deep inelastic scattering (DIS) 
reaction e + p —>■ e + X. The measurements cover the kinematic range Q 2 > 3 GeV 2 and 
45 < W < 225 GeV, where W is the total mass of the photon-proton system. The data for 
x L > 0.93 were used in an earlier study [17] to extract the diffractive structure function 
of the proton. 

The leading proton structure function, F^ p , defined in Section 3, which can be identified 
with a fracture function, is also presented. The latter parameterises the momentum 
spectra of leading particles through parton distribution functions in the proton. This 
approach can be incorporated in Monte Carlo (MC) programs simulating hadronic final 
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states in pp interactions at the LHC [18] and extended cosmic-ray showers [19,20]. 



2 Experimental set-up 

The measurements were performed with data collected in 1997 at the ep collider HERA 
using the ZEUS detector, when HERA operated with a proton beam energy E p = 820 GeV 
and a positron beam energy E e = 27.5 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [21]. A brief outline 
of the components which are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [22], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the polar angle 1 
region 15° < 6 < 164°. The transverse-momentum resolution for full-length tracks was 
°{Pt)/pt = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [23] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = 0.18/V^E for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The position of the scattered positron was determined by combining information from the 
CAL, the small-angle rear tracking detector [24] and the hadron-electron separator [25]. 

The LPS [16] was used during the data-taking period 1994—2000 to detect positively 
charged particles scattered at very small angles and carrying a large fraction of the longi- 
tudinal momentum of the incoming proton. It consisted of 54 planes of silicon microstrip 
detectors grouped into six stations, SI to S6, and located along the outgoing proton beam 
line between Z = 20 m and Z = 90 m. 

The stations SI, S2, S3 and S4, S5, S6 can be considered as two independent spectrom- 
eters, called sl23 and s456, respectively, due to their different phase-space coverage, in 
particular in azimuth. During data taking, the stations were inserted very close to the 
proton beam (typically a few mm). Charged particles inside the beampipe were deflected 
by the magnetic field of the proton-beamline magnets and measured in the LPS with a 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards the 
center of HERA. The coordinate origin is at the nominal interaction point. 
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resolution better than 1% on the longitudinal momentum and of 5 MeV on the trans- 
verse momentum. The beam transverse momentum spread at the interaction point was 
~ 40 MeV in the horizontal plane and ~ 90 MeV in the vertical plane and dominated the 
transverse-momentum resolution. 

A forward neutron calorimeter (FNC) [26] was installed in the HERA tunnel at 9 = 0° 
and at Z = 106 m from the interaction point in the proton-beam direction. The FNC, a 
lead-scintillator calorimeter, had an energy resolution for hadrons a(E)/E = 0.70/v^E, 
with E in GeV, as measured in a test beam. The calorimeter was segmented vertically 
into 14 towers. Three planes of veto counters were located in front of the FNC to reject 
events in which a particle showered in inactive material along the beamline upstream of 
the FNC. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small-angle energetic photons were measured by the luminosity monitor [27], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 



3 Kinematics and cross sections 

Figure 1 illustrates semi-inclusive leading proton production in ep collisions. Four kine- 
matic variables are needed to describe the reaction e + p — > e + Xp. They are defined in 
terms of the four-momenta of the incoming and outgoing positron, k and k', and of the 
incoming and outgoing proton, P and P', respectively. 

The Lorentz-invariant kinematic variables used in inclusive studies are Q 2 = —q 2 = 
— (k — k') 2 , the virtuality of the exchanged photon; x = Q 2 /(2P ■ q) and the inelasticity, 
y = q- P/(k ■ P) ~ Q 2 /(sx); W 2 = (P + k - k 1 ) 2 = m 2 p + Q 2 (l - x)/x, the square of the 
photon-proton centre-of-mass energy. In these equations, m p is the mass of the proton 
and a/s = 300 GeV is the e + p centre-of-mass energy. Among these variables, only two are 
independent. 

Two additional variables are required to describe the leading proton. They are chosen as 
the momentum fraction carried by the outgoing proton 

P'-k 
XL= P~k 

and its transverse momentum with respect to the direction of the incoming proton, px- 
In terms of these variables, the square of the four-momentum transfer from the target 
proton is given by 
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t — {P — P ) — —-^ m p , 



2 

Ill 

XL X L 



where the second term is the minimum kinematically-allowed value of \t\ for a given xl- 
The variable t is the square of the four-momentum of the exchanged particle. 

The differential cross section for inclusive e + p — > e + X scattering, in the Q 2 region of this 
analysis, is written in terms of the proton structure function, F 2 (a;,Q 2 ), as 

d 2 a e+p ^ e +x Ana 2 ( y 2 \ 2 

dxdQ 2 =-^{ 1 - y + j) F ^ Q)il + A) > (1) 

where A is a correction that takes into account the contribution of the longitudinal struc- 
ture function, F L , and of electroweak radiative effects. Similarly, the differential cross 
section for semi-inclusive leading proton (LP) production can be written in terms of the 
leading proton structure function, F^^ix, Q 2 , xl,Pt), as 

where A L p is the analogue of A. 

The structure function F% F ^ (x,Q 2 ,Xl,Pt) corresponds to the proton-to-proton fracture 
function M^ p (x, Q 2 , xl,Pt) [15], i.e. the structure function of a proton probed under the 
condition that the target fragmentation region contains a proton with a given xl and p\. 



4 Reconstruction of the kinematic variables 

In the Q 2 range of this analysis, DIS events are characterised by the presence of a scattered 
positron, mostly in RCAL. The scattered positron was reconstructed using an electron 
finder algorithm based on a neural network [28]. 

The properties of the hadronic final state in the central detector were derived using the 
energy flow objects (EFOs) [29] reconstructed from CAL clusters and CTD tracks by 
combining the CTD and CAL information to optimise the resolution of the reconstructed 
kinematic variables. The EFOs were additionally corrected for energy loss due to inactive 
material in front of the CAL. 

The DIS variables x, y, Q 2 and W were obtained by using a weighting method [30], which 
uses a weighted average of the values determined from the electron [31] and double- 
angle [32] methods. The variable y was also reconstructed using the Jacquet-Blondel 
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method [33], which uses information from the hadronic final state to reconstruct the 
event kinematics, and is denoted by i/jb- 

The momentum p = (px,Py,Pz) of the leading proton candidate was determined using 
the LPS. The variable xl was evaluated as xl = Pz/E p and the squared transverse 
momentum as p\ = p\ + p\ . 

5 Data sample and event selection 

During 1997, the ZEUS detector collected an integrated luminosity of 27.8 pb -1 . However, 
the experimental conditions allowed the operation of the LPS only for an integrated 
luminosity of 12.8 pb -1 . Of this sample, 4.8 pb -1 of data were collected with all the LPS 
stations. In the remaining part, only the spectrometer s456 was used. 

Online, a three-level trigger [34] was used. At the third level, the event variables were 
reconstructed with an accuracy close to that obtained offline. Final detector calibration 
and full-event reconstruction were performed offline. 

Two sets of events were selected [35]: the inclusive DIS sample and the LPS sample. For 
a fraction of the inclusive DIS candidate events the trigger was prescaled, thus reducing 
the effective integrated luminosity of the inclusive DIS sample to 1 pb -1 . The selection 
of the LPS sample was performed with a dedicated LPS trigger. 

The presence of a good positron candidate in the CAL was required in the trigger chain 
used to select the inclusive DIS sample. In addition, the following conditions were applied: 

• l-ZVtxl < 50 cm, where Z vtx is the Z coordinate of the event vertex. This cut is needed 
to remove background due to proton beam-gas interactions and cosmic rays; 

• energy of the scattered positron E' e > 10 GeV. The positron position was required to 
be outside the region close to the rear beampipe hole, where the presence of inactive 
material reduced the precision of the energy measurement; 

• the quantity E — Pz, where the energy E and the longitudinal momentum Pz are 
summed over all the EFOs and the scattered positron, in the range 38 < E — Pz < 65 
GeV, to exclude background from photoproduction, proton beam-gas interactions and 
cosmic rays; 

• 2/jb > 0.03 in order to ensure hadronic activity away from the forward direction. 
The following cuts define the kinematic region: 

• Q 2 > 3 GeV 2 , to select DIS events with large virtuality of the exchanged photon; 

• 45 < W < 225 GeV, to ensure a wide kinematic coverage of the hadronic system. 
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The number of events that passed the inclusive DIS selection cuts was 145447. 

The LPS sample was selected as the DIS sample, but requiring in addition the LPS trigger 
and the following conditions: 

• a reconstructed track in the LPS with p\ < 0.5 GeV 2 and xl > 0.32. To reduce 
the sensitivity of the LPS acceptance due to the uncertainties of the location of the 
beampipe elements, a cut was applied to the minimum distance, A P i pc , between the 
track and the beampipe requiring A pipc > 0.04 (0.25) cm for s456 (sl23). The selection 
of tracks with A pot > 0.02 cm, where A pot is the minimum distance of the track from 
the edge of any LPS detector, ensured that the tracks were well within the active 
regions of the silicon detectors; 

• the sum of the energy and the longitudinal momentum of both the energy deposits in 
the CAL and the particle detected in the LPS, E + P z , was required to be smaller than 
1655 GeV; this cut rejected most of the random overlays of DIS events with protons 
from the beam-halo or from a proton beam-gas collision (see Section 8). 

A total of 73275 events survived the above selection criteria, of which 6008 had a track 
in sl23 and 67267 had a track in s456. 

6 Monte Carlo simulation 

To determine the acceptance of the apparatus, inclusive DIS events with Q 2 > 0.5 GeV 2 
were generated using Djangoh [37], which is interfaced to Heracles [38] for electroweak 
radiative effects. In order to study the migration of events from low Q 2 , a sample of 
photoproduction events with Q 2 < 0.5 GeV 2 was generated with Pythia [39]. In the 
MC samples, the hadronic final state was generated using the Matrix Element Parton 
Shower model (MEPS) [40] for QCD radiation and Jetset [41] for hadronisation. The 
diffractive events in Djangoh were generated using the soft-colour-interaction mechanism 
(SCI) [14]. 

All MC events were passed through the standard simulation of the ZEUS detector, based 
on GEANT 3.13 [42], and trigger and through the same reconstruction and analysis pro- 
grams as the data. The simulation included the geometry of the beampipe apertures, the 
magnetic field along the leading proton trajectory and the proton-beam emittance. 

To obtain a good description of the data, it was necessary to reweight the leading proton 
xl and p\ distributions generated by the MC [35,43]; the fraction of diffractive events 
with respect to the total was also reweighted in bins of xl. In particular, the slope 
of the exponential p\ distribution, ranging from 2.5 to 4.5 GeV -2 , was increased by a 
constant value Ab = 3.4 GeV -2 and the xl spectrum was reweighted to a flat distribution 
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below the diffractive peak. The reweighting parameters were chosen according to previous 
measurements [7]. The reweighting preserved the total MC cross section. 

For the LPS sample, the comparison between the data and the sum of the reweighted MC 
samples (Djangoh and Pythia), for the DIS variables and the LPS specific variables, 
is shown in Figs. 2 and 3. The agreement is generally good. The LPS variable Ap ipc is 
not perfectly reproduced by the reweighted MC, but the selection cut applied is far from 
the region in which the disagreement is observed. 

7 Acceptance 

The acceptance was defined as the ratio of the number of reconstructed events in a bin 
to the number of generated events in that bin. This definition includes the effects of 
the geometrical acceptance of the apparatus, its efficiency and resolution, and the event 
selection efficiency. Figure 4 shows the acceptances of the LPS station combinations sl23 
and s456 as a function of xl and p\. The maximum acceptance is 10% in the region 
0.63 < %l < 0.65, 0.05 < pj, < 0.1 GeV 2 for the spectrometer sl23 and 52% in the region 
0.77 <x L < 0.8, r p\ < 0.05 GeV 2 for the spectrometer s456. 

The analysis bins were chosen according to the LPS acceptance, resolution (see Section 2) 
and available statistics. 

For completeness, also shown in Fig. 4 are the three regions of p\ used for the cross-section 
measurements. 

8 Background studies 

The LPS data sample contains three background contributions, 

• non-baryon contributions; 

• overlay events; 

• misidentified low-Q 2 background. 

The LPS had no particle identification capability, but MC studies indicate that most high- 
xl particles in the LPS are protons. The MC expectations were tested with a subsample 
of LPS-tagged events where a neutron candidate was found in the FNC [36]. The neutron 
candidate was required to have a minimum energy deposit of 50 GeV and the total E + P z 
of the event, including the neutron, was required to be below 1750 GeV. A total of 47 
events were found. For this class of events the track in the LPS is most likely either a n + 
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or a K + . Figure 5a shows the ratio Pfnc of the number of events with a track in the LPS 
and a neutron candidate in the FNC to the number of events with a track in the LPS. 
The reweighted MC describes the data well and therefore can be used to subtract the 
background. The fraction R of events with a positive meson reconstructed in the LPS, 
evaluated with the reweighted MC, is shown in Fig. 5b as a function of Xj_,. The fraction 
R in the MC was found to be independent of p%. It is substantial at low xl and falls 
below 10% above xl ~ 0.6. This contribution was subtracted. 

The E + P z spectrum for the beam-halo events was constructed as a combination of 
generic DIS events and a beam-halo track reconstructed in the LPS in randomly triggered 
events. The E + Pz distribution was normalised to the data for E + Pz > 1685 GeV, which 
mainly contain beam-halo events. The background remaining after the E + Pz < 1655 
GeV cut was negligible for xl < 0.9, and reached (8 ± 3)% for xl > 0.98. The expected 
fraction of overlay events was subtracted. 

The acceptance corrections were calculated using the reweighted MC generated with Q 2 > 
2 GeV 2 . The contribution of events which migrate from the region Q 2 < 2 GeV 2 was found 
to be independent of Xl and p\ and equal to (7.3 ± 0.5)%; it was subtracted. 

9 Systematic studies 

The systematic uncertainties were calculated by varying the cuts and by modifying the 
analysis procedure. The stability of the DIS selection was checked by varying the selection 
cuts, 

• the l-Zvtxl cut was varied by ±10 cm; 

• the cut on the scattered positron energy was varied by ± 2 GeV and the width of the 
fiducial region in the rear calorimeter was varied by 0.5 cm in the X and Y directions; 

• the E — P z cut was changed to 35 < E - P z < 68 GeV and 41 < E - P z < 62 GeV; 

• the yjB cut was raised to 0.04. 

The observed changes in the cross section were below 1% and neglected. The variation in 
the LPS selection of the A pipe threshold by ± 0.03 cm and the A pot threshold by ± 0.01 
led to negligible changes in the cross section (<1%). 

The following checks resulted in non-negligible systematic uncertainties of the cross section 
(the mean value is given in brackets): 

• the reweighting parameter Ab was varied by ±0.9 GeV -2 , compatible with the spread 
of b versus x L (+^); 
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• the subtracted fraction of background from tt + and K + reconstructed in the LPS was 
varied by the statistical uncertainty derived from the LPS-FNC data (see Fig. 5) 

for x L < 0.7); 

• the fraction of overlay events to be subtracted from the data was increased and de- 
creased by its statistical uncertainty (^2 0% for £l > 0.9); 

• the uncertainty on the beam optics was evaluated by varying the transverse momentum 
spread of the proton beam in the MC by ±10% [17] and led to a change of typically 
±1.4%. In addition, the positions of the LPS stations were varied to reflect the actual 
position of the stations during the data-taking period. This was done because, in the 
simulation, the MC assumes only one average position (^8%)- ^ n tne diffractive region 
the uncertainty related to the beam optics increased to ±10%. 

The resulting total systematic uncertainty, obtained by adding in quadrature all the in- 
dividual systematic uncertainties (combining positive and negative contributions sepa- 
rately), excluding an overall normalisation uncertainty of 2% from the luminosity mea- 
surement, is shown in the figures as an error band, that includes both a correlated and 
an uncorrelated component. 

10 Results 

All measurements were performed separately with the sl23 and s456 spectrometers, and 
were then combined in a weighted average, using only statistical uncertainties. This 
procedure was repeated for every systematic check. Some measurements are presented 
as normalised to the inclusive DIS cross section, a- mc , determined from the inclusive DIS 
sample described in Section 5. All cross-section measurements are averaged over a given 
bin and quoted at the mean value of the variable in that bin. The measurements for 
xl > 0.93 are presented here in a different kinematic domain than those previously 
published [17]. 

10.1 Transverse-momentum spectra and slopes 

The double-differential cross-section g? 2 <7lp / 'dx^d'p^. as a function of in bins of xl is 
presented in Fig. 6 and given in Table 1. The results obtained with the sl23 and s456 
spectrometers are shown separately. Within the uncorrelated uncertainties, the two sam- 
ples lead to consistent results. The lines shown in Fig. 6 represent the results of a fit of 
an exponential function A-e~ bp T to the combined cross-section gPolp / 'dxidpip. The band 
shows the statistical uncertainty of the fit. The slopes are presented as a function of xl 
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in Fig. 7 and given in Table 2. The slopes show no strong dependence on ii. The mean 
value of the slopes is (b) = 6.76 ± 0.07(stat.)+o;52(syst.) GeV -2 . The measurements of the 
V 2 T slopes at (Q 2 ) = 5.1 GeV 2 and (Q 2 ) = 30.1 GeV 2 , in the range 45 < W < 225 GeV, 
where xl bins were combined, are presented in Fig. 8 and given in Table 3. Also shown 
are the ZEUS 1994 data [16] in the range Q 2 < 0.02 GeV 2 and 176 < W < 225 GeV and 
the ZEUS 1995 [7] data in the range 0.1 < Q 2 < 0.74 GeV 2 and 85 < W < 258 GeV. The 
p\ slopes are independent of the virtuality of the exchanged photon. 



10.2 Longitudinal momentum spectra 

The cross section as a function of xl has been measured in three bins of p\: < p\ < 0.04, 
0.04 < p\ < 0.15 and 0.15 < p\ < 0.5 GeV 2 . The leading proton production rate, 
1/oinc ■ dcr L p/dxL, for the three p\ ranges is shown in Fig. 9 and listed in Table 4. Due to 
the LPS acceptance, the accessible x L range chang GS clS el function of p\, as seen in the 
figure. The rate as a function of xl is approximately flat up to the diffractive peak, where 
it increases by a factor of about six. This behaviour of the cross section as a function of 
xl is essentially independent of p\. 

Since, as discussed in the previous section, also the p\ slopes are independent of xl, the 
cross section as a function of Xl can be extrapolated to the full xl > 0.32 and p\ < 0.5 
GeV 2 range. The measurement of l/a\ nc ■ daip/d-XL as a function of Xl, extrapolated to 
the fullpf, < 0-5 GeV 2 range is shown in Fig. 10 and given in Table 5. For comparison, the 
ZEUS 1995 data [7] with lower Q 2 are also shown. The two measurements are consistent. 

For p\ < 0.04 GeV 2 , the measurement of l/cr inc • da^p/dxL can also be compared to 
previous measurements in the photoproduction regime (Q 2 < 0.02 GeV 2 ) [7]. The com- 
parison is shown in Fig. 11. Due to the low p\ values, the diffractive peak is not accessible 
(see Fig. 4). The photoproduction data tend to lie systematically below the higher-Q 2 
measurement, though within uncertainties the results are consistent. 



10.3 Ratios of leading proton production to inclusive DIS yields 

The rate of leading proton production, r LP ^(x, Q 2 , Xl), in e + p scattering was determined 
according to 

lp (% n 2 ,_ N LPS (x,Q 2 ,x L )A Dls £ Dls 1 i 
' ^ ' XL > N ^(x,Q 2 ) A lps £ lps Ax l ' [ ) 

where N LPS (x, Q 2 ,Xl) is the number of events corresponding to an integrated luminosity, 
£ LPS , with a proton candidate in the LPS in a given (x,Q 2 ,xl) bin and integrated over 
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< pj> < 0.5 GeV 2 , and N Dls (x, Q 2 ) is the number of DIS events corresponding to an 
integrated luminosity, £ DIS , in that (x, Q 2 ) bin. The acceptance Adis was estimated by 
applying only the DIS selection cuts and Alps is the acceptance of the LPS sample. The 
variable Axl is the size of the xl bin. 

The ratio r LF ^ as a function of xl in bins of x and Q 2 is shown in Fig. 12 and given 
in Table 6. The xl range of the measurement is limited to 0.32 < xl < 0.92, as detailed 
studies of the diffraction region were presented elsewhere [17]. The ratio r LP ^ has also 
been measured in the three ranges of p\ and the values are given in Tables 7, 8 and 9. 
The r LP ( 3 ) values are approximately constant over the kinematic range of this analysis, 
independent of the p\ range. 

The rate of leading proton production as a function of x and Q 2 , integrated over 0.32 < 
x L < 0.92 and p\ < 0.5 GeV 2 , r LP(2 )(x, Q 2 ), is shown in Fig. 13 and given in Table 10. 
The ratio r LF ^ is approximately constant as a function of x and of Q 2 . The mean value 
(r LP ( 2 )) = 0.240 ±0.001(stat.)l^(syst.) means that approximately 24% of inclusive DIS 
events have a leading proton in the range 0.32 < xl < 0.92 with p\ < 0.5 GeV 2 . 

The ratio (r LP ( 2 )) averaged over i as a function of Q 2 is shown in Fig. 14 and given 
in Table 11, in the range 0.32 < x L < 0.92, p\ < 0.5 GeV 2 and 45 < W < 225 GeV. 
A mild increase with Q 2 cannot be excluded. To further investigate the Q 2 dependence, 
the rates integrated over 0.6 < xl < 0.97 can be compared to the equivalent rates for 
(Q 2 ) = 0.29 GeV 2 measured by the ZEUS collaboration previously [7]. The result is shown 
in Fig. 14 and included in Table 11. Assuming that the systematic uncertainties have a 
similar origin, dominated by the LPS acceptance, a steady rise with Q 2 is observed. A 
similar effect was observed in leading neutron production [36], where it was attributed to 
absorption and rescattering effects [44], which disappear when the transverse size of the 
virtual photon decreases with increasing Q 2 . 

10.4 The leading proton structure functions 

The ratio r LP ^ can be expressed as the ratio of F 2 F ^ to F 2 , 



r ^ )(l , g2) = ^L_M), (4) 

where is obtained from F 2 LP ^ by integrating over xl and p\. Therefore, the 

values of F 2 LP ^ can be obtained from the measured r LP ( 2 ) and F 2 . The values of F 2 were 
obtained from the NLO ZEUS-S fit parameterisation of the parton distribution functions 
of the proton [45]. 
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The structure- function F 2 is presented in Fig. 15, plotted as a function of x at fixed 
Q 2 values for 0.32 < x L < 0.92, p\ < 0.5 GeV 2 and given in Table 12. The curves in the 
plot show the F 2 parameterisation scaled by the average value (r LP ( 2 )) = 0.24 and the 
bands represent the one-standard-deviation limits of the NLO ZEUS-S parameterisation. 
A very good description of F 2 is obtained, as expected from the approximate x and 
Q 2 independence of r LP ^ . 

10.5 Comparison to leading neutrons 

The rate of leading proton production for p 2 - <0.04 GeV 2 can be compared to the recent 
ZEUS measurement of leading neutrons [36]. The comparison is shown in Fig. 16. In 
the range 0.32 < xl < 0.92, there are approximately twice as many protons as neutrons. 
This is consistent with the additive quark model [46], in which the probabilities to have 
a proton or a neutron in the final state are 2/3 and 1/3, respectively. In a particle 
exchange model, the exchange of isovector particles would result in half as many protons 
as neutrons. Thus, exchange of isoscalars must be invoked to account for the observed 
proton rate [13]. This contribution is likely also to explain the different behaviour of the 
rates at large xl. 

The slopes of the p\ distributions for leading protons and neutrons are shown in Fig. 17. 
Although the p\ and Q 2 ranges are different, the two samples have similar values of b in 
the region 0.65 < xl < 0.8, where pion exchange is expected to dominate the production 
mechanism [13]. 

10.6 Comparison to models 

The predictions of the model of Szczurek et al. [13] are compared to the leading proton 
production rate 1/ a mc -da\jp / dxi and the p\ slopes in Fig. 18. In this model, leading proton 
production for 0.6 < Xl < 0.9 is dominated by isoscalar Reggeon exchange. Diffractive 
processes due to Pomeron exchange become increasingly important as il approaches 
unity. The contribution of pion exchange plays an important role in the medium xl 
range. The model describes the shape of the longitudinal momentum spectrum and of 
the pfn slopes reasonably well. The normalisation of the Reggeon contribution, which 
has a large theoretical uncertainty [13], may be constrained by this measurement. The 
model does not include absorptive corrections and rescattering effects [44], since they are 
expected to be small in DIS regime. 

In Fig. 19, various DIS Monte Carlo models are compared to the experimental data. 
The prediction of Djangoh [37] with SCI [40] and Rapgap [47] are compared to the 
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leading proton production rate l/c inc • da^/dxi and to the p\ slopes. In both MC 
models, the QCD radiation was performed either by the parton shower [40] or colour 
dipole (CDM) [48] models. None of the DIS Monte Carlo models can reproduce the flat 
dependence of xl below the diffractive peak. The MC generator Djangoh, with SCI and 
MEPS, reproduces quite well the dependence of b on Xl, although the mean values of the 
slope are lower than those measured. In the other MC models, the value of the slope is 
consistent with the measurements only at high xl. 

11 Summary 

The cross section of leading proton production for xl > 0.32 and p\ < 0.5 GeV 2 and its 
ratio to the inclusive DIS cross section have been measured in the range Q 2 > 3 GeV 2 
and 45 < W < 225 GeV, with 12.8 pb" 1 collected during 1997 with the ZEUS leading 
proton spectrometer. The leading proton production cross section as a function of p\ falls 
exponentially with a mean slope (b) = 6.76±0.07(stat.)^Q;52(syst.) GeV -2 , approximately 
independent of xl and of the photon virtuality, Q 2 . Below the diffractive peak, the xl 
distribution is flat, independent of p\ and Q 2 . 

The ratios of leading proton production to the inclusive DIS yields, r LP ^ and r LP ( 3 \ show 
no strong dependence on x or Q 2 . In the range 0.32 < xl < 0.92 and p\ < 0.5 GeV 2 , 
approximately 24% of DIS events have a leading proton. The ratio (r LP ( 2 )) averaged 
over x rises very slowly with Q 2 in the DIS regime. This trend is further confirmed 
by measurements at lower Q 2 . The dependence of the leading proton structure- function 
F 2 LP( - 2 ' ) on x and Q 2 is approximately the same as that of F 2 . 

The yield of leading protons in DIS is almost twice that of leading neutrons. The p\ slopes 
have a different dependence on xl, but are compatible in the range 0.65 < xl < 0.8, where 
pion exchange is dominant. 

The main features of the experimental data are reproduced by a Regge-inspired model. 
The results of this study are an important ingredient for modelling leading particle pro- 
duction in hadron-hadron interactions, which is not properly reproduced by existing gen- 
erators. 
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0.15-0.20 




136 ± 22^n 


136 ± 22±^ 


0.62-0.65 


0.00-0.05 


410 ± 53+^ 


463 ± 20±i? 


456 ± 19±\l 




0.05-0.10 


1 1 Q 

296 ± 29±i| 


i 1 ft 

324 ± 3llgf 


I 1 O 

309 ± 21±H 




0.10-0.15 


189 ± 43tH 


187 ± 16±\l 


187 ± 15±\ 2 S 




0.15-0.20 




148 ± 17if 7 


148 ± 17+Jt 




0.20-0.25 




134 ± 22+^ 2 


134 ± 22+^ 2 


0.65-0.68 


0.00-0.05 


421 ± 74+^ 


464 ± 16±?| 


462 ± 16±H 




0.05-0.10 


320 ± 29±H 


297 ± 33±§| 


310 ± 22±H 




0.10-0.15 


181 ± 36^g 


298 ± 37l^ x 


238 ± 26128 




0.15-0.20 




118 ± 13lg 2 


118 ± 13lf 




0.20-0.25 




119 ± 16lg 8 


119 ± 16^g 8 




0.25-0.35 




87 ± 12±X 5 


87 ± 12± r 15 



Table 1: The double- differential cross-section d 2 o~up / dxLdp% as a function of xl 
andp\, separately measured with the spectrometers sl23 and sJf.56 and the combined 
result. Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Table 1 (cont.) 



xl range 


p\ range 
(GeV 2 ) 


d 2 crup 1 'dx Ldpr 
si 23 (nb/GeV 2 ) 


d 2 a\jp / dx idpip 
s456 (nb/GeV 2 ) 


d 2 (T\jp 1 'dx^dprn 
Combined (nb/GeV 2 ) 


0.68-0.71 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


253 i 851^34 
320 ± 29^42 
184 ± 31±f^ 
92 ± 41^° 


500 ± 1S±\1 
366 ±41^ 8 
259 ± 47±|§ 
155 ± 24^? 
103 ± 16lf 
85 ± 10±2 

44 ± 9tg 


493 ± 15±H 
335 ± 24±\l 
207 ± 26±^ 
139 ± 21±H 
103 ± 16lf 
85 ± 10±| 

44 ± 9lg 


0.71-0.74 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


279 ± 27±f§ 
217 ± 33iy 
114 ± 43±H 


472 ± 12^11 
255 ± 24+1^ 
130 ± 33±?| 
250 ± 661^3 
116 ± 27+f| 
86 ± 16±H 
27 ± 51 2 . 


472 ± 12+^| 
266 ± 181^ 
173 ± 23±f§ 
154 ± 36^69 
62 ± 19±f 5 
86 ± 16±^ 
27 ± 5±l 


0.74-0.77 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


250 ± 30±|^ 
234 ± 30±| 3 
108 ± 35±|f 


467 ± 10t^ 8 
325 ± 21^ 2 
199 ± 40±H 

152 ± 52±|| 
121±34l|o 
8 ± 41} 1 


467 ± lO^g 
300 ± 17tfl 
222 ± 24lf 
108 ± 35±|f 
152 ± 52±|| 
121 ±341°, 
8 ± A+l 1 


0.77-0.80 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 


248 ± 39tj$ 
196 ± 23lf 
198 ± 48±|§ 


452 ± 10±$ 
297 ± 14±^| 
241 ± 20^1° 
203 ± 32± 9 9 7 
168 ± 54±|| 


452 ± 10l}? 
292 ± 13l^ 
222 ± 15±| 2 
201 ± 27±f| 
168 ± 54±|^ 


0.80-0.83 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 


353 ± 74l JJ 7 
197 ± 25+f2 
84 ± 23l 2 ^ 
73 ± 32 j|| 


451 ± 10l? 9 
266 ± lllg 8 
210 ± Ut 7 16 
178 ± 17^ 5 
136 ± 16±§ 3 


451 ± 10l? g 
268 ± lllg 8 
207 ± 12t^ 3 
147 ± 14±jj 3 
123 ± 15±1| 
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Table 1 (cont.) 



xl range 


p\ range 
(GeV 2 ) 


d 2 (j\jp / dx Lctptp 
sl23 (nb/GeV 2 ) 


d 2 (Tip / dx idpip 
s456 (nb/GeV 2 ) 


d 2 ; clp I 'dx hdp\ 
Combined (nb/GeV 2 ) 




0.25-0.35 




98 ± 14±f 


98 ± 14± 22 


0.83-0.86 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


184 ± 24+Jo 
190 ± 37+jg 5 
50 ± 20+J5 


434 ± 10±i| 
295 ± 12±§3 
207 ± 13±f § 
169 ± 15±| 9 
127 ± 14lf 
122 ±14+^ 
49 ± 8lg 2 


434 ± 10±i| 
295 ± 12+J 3 
202 ± 11+^ 7 
172 ± 14±| 9 
102 ± H±i9 
122±14l 3 7 6 
49 ± 8± 6 2 


0.86-0.89 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


247 ± 36±H 
196 ± 32^ 
115 ±34+^ 
24 ± 16+^ 


467 ± 13+f| 
309 ± 12±|g 
232 ± 14+J 
180 ± 15i^ 5 
136 ± 14+^ 
81 ± 8+55 
43 ± 6+J 


467 ± 13+5g 
309 ± 12±|g 
234 ± 13±J 9 
183 ± 14+^ 
133 ± 13+5° 
69 ± 7+g 4 
43 ± 6+J 


0.89-0.92 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


262 ± 47lff 
172 ± 27±|| 
97 ± 27+^s 
37 ± 19+55 


481 ± 19±f| 
315 ± 13±^ 
228 ± U±H 
175 ± U+H 
133 ± 13t\l 
69 ± 7t\ 
39 ± 5+J 


481 ± 19±f| 
315 ± 13±^g 
231 ± 13la4 
174 ± 12^55 
126 ± 12+55 
65 ± 6+J 
39 ± 5+J 


0.92-0.95 


0.00-0.05 
0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 
0.35-0.50 


221 ± 58±§J 2 
138 ± 21+^ 
92 ± 22t5o 
39 ± 15+5*8 


820 ± U5±\ji 
235 ± 14±| 4 
182 ± UtH 
156 ± 14+5g 
114 ± 13+io 
83 ± 9+5^ 
36 ± 5tf 


820 ± 145+^ 
235 ± 14±| 4 
184 ± 13±|f 
150 ± 12+J° 
109 ± ll±if 
72 ± 8+5 3 
36 ± 5±f 


0.95-0.98 


0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 


280 ± 39l5e 
254 ± 46l5i6 
50±17±Jg 


423 ± 35+^6 
328 ± 30l5o 
250 ± 32+J 42 
197 ± 32+J° 
190 ± 34t?2 


423 ± 35+J? 6 
328 ± 30l5o 
262 ± 25+5° 
215 ± 26^26 
77 ± 15±i» 
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Table 1 (cont.) 



xl range 


p\ range 
(GeV 2 ) 


d 2 (j\jp / dx Lctptp 
si 23 (nb/GeV 2 ) 


d 2 a\jp / dx idpip 
s456 (nb/GeV 2 ) 


d 2 a\jp 1 'dx^dplj, 
Combined (nb/GeV 2 ) 


0.98-1.00 


0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.35 


1218 ± 1601^ 
574 ± 8lt^ 6 
231 ± 54±^ 6 


2788 ± 180±f^| 
1423 ± 90±i|f 
1012 ± 89±i|§ 
848 ± 93l|f 
588 ± 69ti 


2788 ± 180^ 
1423 ± 90t\ll 
1061 ± 78±ff 8 
693 ± 6ll^ 4 
367 ± 42l^ 7 
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Xl ranj 


? e 


o (GeV ) 





50 - 





56 


U - JO -1.06-1.53 





56 - 





62 


8 4q+0.47+2.18 
°-^ J -0.51-1.32 





62 - 





65 


. i_n ia -L 1 n 

7 q/?H-U.o4H-l.lz 

' - OU -0.38-0.54 





65 - 





68 


_ i_n or_l n 7S 

f} oqH-U.zoH-U. f o 

U.OO_o,26-0.38 





68 - 





71 


, , in 01 i n i k 

7 m +u.^i+u.io 

' - ui -0.21-0.51 





71 - 





74 


' -^°-0.26-0.38 





74 - 





77 


7 A Q+0.OO + 1.ZO 

' -^°-0.38-1.00 





77- 





OA 

oU 


fi OO+0.30+0.90 
U - OZ -0.34-1.24 





80 - 





83 


e c^+0.21+0.52 
u - d ^-0.21-0.40 





83 - 





86 


c OO+0.17+0.75 
J - oo -0.17-0.18 





86 - 





89 


fi OO+0.17+0.82 
U - OZ -0.17-0.64 





89 - 





92 


7 99+O.2I+O.72 
' - zz -0.21-0.75 





92 - 





95 


fi nn +0.26+0.48 
u,uu -0.26-1.37 





95 - 





98 


a 4 r+0.47+1.69 
^•^ d -0.47-0.88 





98 - 


1 


00 


o qi +0.34+1.05 
°- oi -0.38-0.91 



Table 2: The p\- slope, b, of the cross-section gPctlp / ' dx^dp\, as defined by the pa- 
rameterisation A-e~ b ' PT and obtained from a fit to the data in bins of Xl. Statistical 
uncertainties are listed first, followed by systematic uncertainties. 



Xl range 


b (GeV" 2 ) 




(Q 2 ) = 5.1 GeV 2 


(Q 2 ) = 30.6 GeV 2 


0.5 - 0.65 
0.65 - 0.8 
0.8 - 0.92 
0.92 - 1.0 


« qn+o.38+i.34 

u - JU -0.43-0.18 

7 nn+°- 17 + - 20 

' - uu -0. 17-0.21 
7 9Q+0.13+0.15 
' - zy -0.17-0.49 
a 7A +0.26+0.28 
°- ' ^-0.30-0.33 


7 Q7+0.43+1.06 
' - J ' -0.38-0.39 
7 1 q+0.17+0.16 
' - io -0.17-0.31 
c i n +0.13+0.64 
u - lu -0.13-0.04 
c nfi +0.30+0.66 
u - uu -0.30-0.36 



Table 3: The p%- slope, b, of the cross-section d 2 o~\xjdxLdp\, as defined by the 
parameterisation A-e~ b ' PT and obtained from a fit to the data in bins ofxi, measured 
in two ranges of Q 2 . Statistical uncertainties are listed first, followed by systematic 
uncertainties. 
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x l ranj 


? e 










l/0"inc 


• da LP /dx L 












< p\ < 0.04 GeV 2 


0.04 < p% 


C 0.15 GeV 2 


0.15 < 


p% < 0.5 GeV 2 


0.32 - 





38 





091 ±0 


nn8 +0.009 
UU °-0.013 
















0.38 - 





44 





082 ±0 


nf)4 +0.005 





156 ± 


016+ 0013 

U.U1U_ 024 










0.44 - 





50 





086 ±0 


nn4 +0.003 





144 ± 


n nio+0.005 
U.U1-3_ 0021 










0.50 - 





56 





091 ±0 


nf)4 +0.004 





140 ± 


0.006±8^g 










0.56 - 





62 





097 ±0 


nn4 +0.002 





143 ± 


0.0041H?! 










0.62 - 





68 





096 ±0 


nn o+0.003 
uuo -0.005 





142 ± 


005+ 002 
u.uuo_ 021 





126 


±0 


nn7 +o.oio 

uu ' -0.011 


0.68 - 





74 





100 ±0 


nn9 +0.003 
uuz -0.005 





141 ± 


005+ 004 
u.uud_ 012 





138 


±0 


nn o+0.008 
UU °-0.009 


0.74 - 





80 





094 ±0 


nn i+o.ooo 

uui -0.007 





145 ± 


nns +0 - 003 





155 


±0 


ni cr+0.010 
uiJ -0.015 


0.80 - 





86 





091 ±0 


nni +o.ooi 

uui -0.005 





139 ± 


03+ - 001 
u.uuo_ 0008 





146 


±0 


QQC-+0.010 

UUJ -0.023 


0.86 - 





92 





099 ±0 


nn9 +0.002 
uuz -0.007 





154 ± 


nns+ - 005 





143 


±0 


nn4 +o.on 

UU4 -0.014 


0.92 - 





98 











167 ± 


006+ - 020 
u.uuu_ 015 





179 


±0 


nf)7 +0.021 
uu ' -0.013 


0.97- 


1 


00 








1 


126 ± 


_o.092 





816 


±0 


nQ«+0.088 
UOU -0.055 



Table 4: The leading proton production rate, l/Vinc • da^/dxL, as a function of 
xl measured in three ranges of p\. Statistical uncertainties are listed first, followed 
by systematic uncertainties. 
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Xl ran; 


?e 


l/0"inc • 


da LP /dx L 





32 - 





38 





372 


± 





n9fi +0.046 
uzu -0.062 





38 - 





44 





309 


± 





n i O+0.037 
ulo -0.037 





44 - 





50 





339 


± 





nl 9+0.033 
ulz, -0.044 





50 - 





56 





358 


± 





nin+0.047 
uiu -0.030 





56 - 





62 





371 


± 





nnq +0.029 
uuy -0.039 





62 - 





65 





371 


± 





m i +0.025 
uii -0.037 





65 - 





68 





385 


± 





nln +0.027 
Uiu -0.031 





68 - 





71 





418 


± 





nin+0.030 
uiu -0.037 





71 - 





74 





398 


± 





nnq+o.039 

UUJ -0.036 





74 - 





77 





408 


± 





nn o+0.038 
UU8 -0.039 





77 - 





80 





398 


± 





nn7 +0.036 
uu ' -0.037 





80 - 





83 





376 


± 





006+ - 029 

uuu -0.035 





83 - 





86 





382 


± 





nnfi +0.026 
UUO -0.033 





86 - 





89 





405 


± 





nn? +0.024 
uu ' -0.030 





89 - 





92 





395 


± 





nna +0.019 
UU8 -0.026 





92 - 





95 





325 


± 





qiq+0.044 
Uiu -0.016 





95 - 





98 





562 


± 





noq+0.058 
UZO -0.049 





98 - 


1 


00 


2 


478 


± 





n 7fi +0.235 
u ' u -0.136 



Table 5: The leading proton production rate, l/a mc ■ da^/dxi, as a function of 
Xl measured in the region p\ < 0.5 GeV 2 . Statistical uncertainties are listed first, 
followed by systematic uncertainties. 
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( 


x) 




(<9 ) (GeV 2 ) 


xl range 




9.6 • 


10" 


•5 


4.2 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.72 
0.82 - 0.72 


0.365 ± 0.086lJu68 
0.268 ± 0.040±g;g|5 
0.391 ± 0.04lio;o4B 
0.338 ± 0.026±g;g|| 
0.378 ± 0.024±g;g|g 
0.458 ± 0.027lg;g^ 


1.7- 


10" 


-4 


4.2 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.414 ±0.063ig;53° 
0.310 ± 0.032±g;g|| 
0.370 ± 0.026lg;g^ 
0.366 ± 0.019±g;g§| 
0.390 ± 0.0151°;°^ 
0.378 ± 0.014±°;°! 


3.5 • 


10" 


-4 


4.2 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.351 ± 0.061±g;g|| 
0.269 ± 0.029±g;g|| 
0.354 ± 0.025lg;gg 
0.406 ± 0.022±°;°|° 

0.405 ± 0.017lg;o38 
0.364 ± 0.014±g;g|| 


6.9 • 


10" 


-4 


4.2 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.275 ± 0.057±g;g|| 
0.308 ± 0.035±g;g|f 
0.339 ± 0.027lg;o34 
0.326 ± 0.020±g;g|| 
0.400 ± 0.018^;^ 
0.403 ± 0.017±g;gfg 


1.46 


• 10 


-3 


4.2 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.272 ± 0.083^;^ 
0.268 ± 0.045;°;°! 
0.285 ± 0.038lg;g^ 
0.355 ± 0.0291^3 

0.362 ± 0.022±g;g|2 
0.359 ± 0.02lt°;°^ 



Table 6: The leading proton production rate, r LP ( 3 \ measured as a function oJxl 
for protons with p\ < 0.5 GeV 2 , in bins of x and Q 2 , with averages (x) and (Q 2 ) . 
Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Table 6 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


1.9 • 10" 4 


7.3 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.314 ± 0.084+°;°^ 
0.430 ± 0.061±g;?|| 
0.299 ± 0.031^^2? 
0.355 ± 0.027±^9 
0.382 ± 0.024to;ol2 
0.383 ± 0.022±g;g|| 


3.4 • 10" 4 


7.3 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.340 ± 0.076±g;g^ 
0.360 ± 0.044ig;g^ 
0.369 ± 0.032±°;° 44 
0.347 ± 0.022±° ) ;° 1 3 ° 
0.374 ± 0.019lg;g^ 
0.359 ± 0.016lg;g^ 


6.9 • 10~ 4 


7.3 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.392 ± 0.0811°;°^ 
0.261 ± 0.038i°;°|? 

0.372 ± 0.034±g:Sai 
0.363 ± 0.025±g:g|| 
0.385 ± 0.020lg;g^ 
0.388 ± 0.019lg;o3? 


1.36 • 1(T 3 


7.3 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.275 ± 0.062io:oei 
0.401 ± 0.056±g;g|g 
0.393 ± 0.039±g;g|| 
0.448 ± 0.0321°;°^ 
0.394 ± 0.021±g;{$§ 
0.389 ± 0.020j£o25 


2.67 • 10~ 3 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.290 ± 0.097lg;^ 
0.442 ± 0.086±g;?§| 
0.449 ± 0.050+^;^6 
0.398 ± 0.033±g;g|| 


2.6 • 10~ 4 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.349 ± 0.078j£?28 
0.285 ± 0.034+°;^ 
0.485 ± 0.037t^74 
0.371 ± 0.021 to^l 
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Table 6 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


4.6 • 10" 4 


11 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.482 ± 0.09718;?^ 
0.285 ± 0.036lg;os9 
0.370 ± 0.031±J];^ 4 
0.435 ± 0.028^054 
0.386 ± 0.019±g;g|J 
0.365 ± 0.016±g;gf| 


9.2 • 10" 4 


11 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.418 ± 0.084±g:g|| 
0.286 ± 0.037±g;g|| 
0.408 ± 0.034lao35 
0.392 ± 0.024l°;°p 
0.392 ± 0.019±g;g|§ 
0.397 ± 0.018±g;g|| 


1.83 • 1(T 3 


11 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.361 ± 0.075lg:o56 
0.353 ± 0.047tjj;i)77 
0.375 ± 0.033±g;g?| 
0.411 ± 0.026±g;g|? 
0.348 ± 0.017±g;g|| 
0.378 ± 0.017±g;g|i 


3.98 • 1(T 3 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.228 ± 0.063±g;g|g 
0.378 ± 0.047lg;g^ 
0.421 ± 0.033±g:2o5 
0.448 ± 0.026lg;oi? 


5.1 • 10" 4 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.470 ±0.11llg;g|J 
0.306 ± 0.037±g;g|| 
0.436 ± 0.035;^4 4 
0.388 ± 0.023lg;SS 


9.2 • 10" 4 


22 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.356 ± 0.084±g;J?J 
0.350 ± 0.053±g;g^| 
0.376 ± 0.039±g;g|g 
0.401 ± 0.030lg;^7 
0.386 ± 0.022lg;gy 
0.411 ± 0.0221°;°^ 


1.84 • 10~ 3 


22 


0.32 - 0.42 
0.42 - 0.52 


0.367 ± 0.0841°;°^ 
0.241 ± 0.038lg;g^ 
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Table 6 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 






0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.347 ± 0.034!°,;°,! 
0.353 ± 0.025±g;gf? 

0.378 ± 0.021±g:o35 
0.401 ± 0.0211°;°^ 


3.66 • 1(T 3 


22 


0.32 - 0.42 
0.42 - 0.52 
0.52 - 0.62 
0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.297 ± 0.078±°,;oi 
0.386 ± 0.066t°;o34 
0.457 ± 0.055±g:gff 
0.442 ± 0.038±g;g^ 
0.453 ± 0.029±g;g|§ 
0.377 ± 0.0231°;°^ 


7.83 • 10~ 3 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.215 ± 0.066l8;^2 
0.354 ± 0.053±g;g|J 
0.399 ± 0.0361°;°^ 
0.442 ± 0.0301°;°^° 


1.03 • 10~ 3 


44 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.356 ±0.104lg;J^ 
0.327 ± 0.0521°;°^ 
0.373 ± 0.038l{];^ 
0.381 ± 0.030±g;g§5 


1.86 • 10~ 3 


44 


0.37 - 0.42 
0.47 - 0.52 
0.57 - 0.62 
0.67 - 0.72 
0.77 - 0.82 
0.87 - 0.92 


0.508 ± 0.150±g;i?e 
0.346 ± 0.073±g;g^ 
0.444 ± 0.063l{];?^ 
0.348 ± 0.0411°;°^ 

0.413 ± 0.0331^52 
0.352 ± 0.026l°,;°J} 


3.68 • 10~ 3 


44 


0.37 - 0.42 
0.47 - 0.52 
0.57 - 0.62 
0.67 - 0.72 
0.77 - 0.82 
0.87 - 0.92 


0.221 ± 0.088±g;g|| 
0.456 ± 0.097±g;22i 
0.303 ± 0.046±g;g|| 
0.496 ± 0.057lHIo 
0.511 ± 0.04418:8" 
0.408 ± 0.03318:829 


7.33 • 10~ 3 


44 


0.37 - 0.42 
0.47 - 0.52 
0.57 - 0.62 
0.67 - 0.72 


0.628 ± 0.226l8;i8o 
0.341 ± 0.090l8;8eo 
0.418 ± 0.074l°;J| 
0.458 ± 0.06118:870 
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Table 6 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 






0.77 - 0.82 
0.87 - 0.92 


0.501 ± 0.048±g;g|f 
0.448 ± 0.039lg;gf§ 


1.54 • i(r 2 


44 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.453 ±0.1351°;°™ 
0.454 ± 0.062±g;g^f 
0.385 ± 0.033±°,;°|f 


2.00 • 10~ 3 


88 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.328 ±0.162+°;^ 
0.462 ±0.12llg;J^ 
0.307 ± 0.045lg;^ 
0.339 ± 0.040±g;g|| 


3.59 • 10~ 3 


88 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.449 ± 0.1451^34 
0.509 ± 0.090+°;J 2 ° 
0.402 ± 0.049±g;J§| 
0.393 ± 0.034±g;g|| 


7.37 • 10~ 3 


88 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.491 ± 0.1781°;^ 
0.345 ± 0.078lg;g^ 
0.364 ± 0.046±g;g§| 
0.440 ± 0.042l{j;of2 


1.42 • 10" 2 


88 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.351 ± 0.133t[];o54 
0.479 ± 0.101±g;}g§ 
0.390 ± 0.054±g;g§f 
0.498 ± 0.050lg;m 


3.01 • 10~ 2 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.449 ± 0.211+°;°^ 
0.330 ± 0.063;2;o 2 3 
0.440 ± 0.056lg;gt| 


4.00 • 10~ 3 


237 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.250 ±0.153lg;u9 
0.175 ± 0.068lg:m3 

i n nvo 

0.451 ± 0.07llg;o44 


7.52 • 10~ 3 


237 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.467 ±0.1951^ 
0.482 ± 0.088lg;^ 
0.430 ± 0.049±g;g|§ 


1.47- 10~ 2 


237 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.301 ± 0.156±g:o40 
0.249 ± 0.074+°;°^ 
0.445 ± 0.070+°;J 4 2 
0.381 ± 0.0451°,;°,^ 
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Table 6 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


3.25 • 10" 2 


237 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.102 ± 0.082±{j;£}? 
0.355 ± 0.090±g;ora 
0.359 ± 0.055tg;o28 
0.403 ± 0.049lg;gg 
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(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


9.6 • icr 5 


4.2 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.082 ± 0.015+5];^ 
0.091 ± 0.012lg;gJ| 

0.080 ± 0.006±o:o26 
0.104 ± 0.008+^;^ 


1.7- 10" 4 


4.2 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.103 ± 0.013^;^ 
0.083 ± 0.008ig;g?g 
0.096 ± 0.005±g;g§| 
0.089 ± 0.004±g:gg| 


3.5 • 10" 4 


4.2 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.077 ±0.010+^ 

0.087 ± 0.008±g;gS§ 

0.108 ±0.0O6lg;oo? 
0.093 ± 0.004±g;gg| 


6.9 • 10" 4 


4.2 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.077 ±0.012+°;°^ 
0.085 ± 0.009lg;g^ 
0.095 ± 0.006±g;gg| 
0.087 ± 0.004lg;g^ 


1.46 • 10" 3 


4.2 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.074 ± 0.016±g;8i§ 
0.046 ± 0.008±g;g^ 
0.083 ± 0.006^^03 
0.086 ± 0.006±g:gg| 



Table 7: The leading proton production rate, r LP ^ , measured as a function of 
xl for protons with p\ < 0.04 GeV 2 , in bins of x and Q 2 , with averages (x) and 
(Q 2 ) . Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Table 7 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


1.9 • 10" 4 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.093 ± 0.017+°;°°;? 
0.098 ± 0.012±°;°°s 
0.087 ± 0.007lg;{j£jj 
0.092 ± 0.007±g;g2i 


3.4 • 10" 4 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.090 ± 0.015lg;gJ| 
0.103 ± 0.011±g;g5S 
0.079 ± 0.005±g;g^ 
0.088 ± 0.005±E];^| 


6.9 • 10" 4 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.090 ± 0.015±g;^| 
0.093 ± 0.011±g;g?f 
0.102 ±0.0071^6 
0.095 ± 0.006±g;g^ 


1.36 • 1(T 3 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.075 ±0.0141°;°^ 
0.104 ±0.013±g;gJ2 
0.113 ± 0.008j£oi7 
0.098 ± 0.006±g:gg| 


2.67 • 10" 3 


7.3 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.067 ± 0.025lg:ral 
0.112 ±0.028lg;g^ 
0.097 ± 0.015±ffi 
0.095 ± 0.01ll^?| 


2.6 • 1(T 4 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.098 ± 0.024t°;°°^ 
0.075 ± 0.013±g:gg| 
0.129 ±0.0141°;°°] 
0.085 ± 0.008±g;gg| 


4.6 • 10" 4 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.114 ±0.018±2;o27 
0.081 ± 0.0091^^4 
0.101 ± 0.007±g;g?f 
0.099 ± 0.006±g:gg| 


9.2 • 1(T 4 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.080 ± 0.013±g:ooi 
0.110 ±0.0121^ 
0.080 ± 0.006±g;gg? 
0.091 ± 0.005lg;g^ 
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Table 7 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


1.83 • 10" 3 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.082 ± 0.0151q;oo7 
0.094 ±0.01llg;g?f 
0.100 ±0.007j£oo6 
0.089 ± 0.005j£ijo2 


3.98 • lO" 3 


11 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.041 ± 0.015^020 
0.118 ±0.021+°;^° 

0.102 ±0.010lK 
0.107 ± O.0O9lg;g54 


5.1 • 10" 4 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.098 ± 0.029±g;g|| 
0.081 ± 0.014±g:gjg 
0.114 ±0.014+°;°^ 
0.088 ± 0.009±g:g§| 


9.2 • 10" 4 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.094 ± 0.018lg;g^o 
0.094 ± 0.012±g;gJ| 
0.100 ± 0.008j£oo8 
0.090 ± 0.006lg:oo4 


1.84 • 10~ 3 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.089 ± 0.017±g;gg| 

0.080 ± 0.010±g;g?5 
0.089 ± 0.007lg;g{$ 
0.089 ± 0.006lg:oo4 


3.66 • 10~ 3 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.078 ± 0.017lg:So5 
0.120 ±0.019lK 

0.114 ± O.OlOloioio 
0.099 ± 0.007l^oi 


7.83 • 10~ 3 


22 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.054 ± 0.019±g;g|2 
0.087 ± 0.018±g;g^ 

i n r\A 1 

0.087 ± o.oioi[j;[] 4 5 

0.116 ± 0.012±g;gf| 


1.03 • 10~ 3 


44 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.092 ± 0.03llg;oi? 
0.086 ± 0.020lg;g^ 
0.105 ±0.016lo:o3i 
0.108 ± 0.013lg;oo5 


1.86 • 10" 3 


44 


0.32 - 0.47 
0.47 - 0.62 


0.096 ± 0.023+°;^ 4 
0.087 ± 0.016±g:g§| 
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Table 7 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 






0.62 - 0.77 
0.77 - 0.92 


0.130 ± 0.015+°;°^ 
0.102 ± 0.010lg;g?| 


3.68 • 1(T 3 


44 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.077 ± 0.023l{j:oii 
0.079 ± 0.016ig;g^ 
0.139 ±0.017±8;|)23 
0.117 ±0.012t°;°° 3 


7.33 • 10~ 3 


44 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.124 ± 0.039±Jj;[!3o 
0.104 ±0.023ig;gg£ 

0.107 ± 0.015tg:oo4 
0.096 ± 0.011±g;gg| 


1.54- 1(T 2 


44 


0.32 - 0.47 
0.47 - 0.62 
0.62 - 0.77 
0.77 - 0.92 


0.076 ± 0.036+°;°^ 
0.118 ±0.034±°;° 2 ° 
0.108 ± 0.019±g;oo9 
0.101 ± 0.015lgioM 


2.00 • 10~ 3 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.130 ± 0.053±g;g^e 
0.128 ± 0.031±g;g|| 
0.089 ± 0.014t°;° 2 ^ 


3.59 • 10~ 3 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.126 ± 0.037±g;g2i 
0.080 ± 0.014t°;° 3 ^ 

0.104 ± 0.012lg;ooi 


7.37 • 10~ 3 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.091 ± 0.031+^;^5 
0.077 ± 0.017+^ 
0.090 ± O.Olllg;^ 


1.42 • 10" 2 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.075 ± 0.024±g;g^| 
0.118 ±0.0271^ 
0.103 ± 0.014t°;° 3 ^ 


3.01 • 10~ 2 


88 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.151 ± 0.075tjj;? 22 
0.081 ± 0.020lg;g^ 
0.099 ± 0.017±g;gg| 


4.00 • 10~ 3 


237 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.053 ± 0.036lg:M5 
0.092 ± 0.033lo;Jo6 
0.116 ±0.025t^ 2 o 


7.52 • 10~ 3 


237 


0.32 - 0.52 
0.52 - 0.72 


0.116 ± 0.052l°;° 3 ^ 
0.108 ± 0.030±g;gg? 
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Table 7 (cont.) 



(x) 


(O 2 ) (GeV 2 ) 




r LP(3) 






0.72 - 0.92 


0.098 ± 0.017l°;™J 


i.47- i(r 2 


237 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.072 ± 0.036ig;g^ 
0.098 ± 0.025±g;g|| 
0.128 ± 0.020j£jj2§ 


3.25 • 10- 2 


237 


0.32 - 0.52 
0.52 - 0.72 
0.72 - 0.92 


0.033 ± 0.020±g;g5i 
0.071 ± 0.019±g;gf| 
0.102 ± 0.017±g:oi9 
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(x) 


(Or) (GeV 2 ) 


xl range 




9.6 • 10~ 5 


4.2 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.146 ± 0.027j£o3g 
0.125 ±0.012t°;°Jf 
0.177 ±0.013t°;°£ 


1.7- 10" 4 


4.2 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.133 ± 0.017^26 
0.145 ± 0.0101K 
0.137 ±0.0O6l8;j)O5 


3.5 • 10" 4 


4.2 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.123 ± 0.018i£oi? 
0.140 ±0.009lg;g?§ 
0.141 ± 0.007±g:ow 


6.9 • 10" 4 


4.2 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.128 ± 0.019l8;o28 
0.128 ±o.oiol8;™8 
0.151 ± 0.008tJj:o2o 


1.46 • 10" 3 


4.2 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.160 ± 0.033l8;ofe 
0.156 ± 0.018±g;gi| 
0.136 ± 0.0101^20 


1.9 • 10" 4 


7.3 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.126 ± 0.023l8;oo6 
0.123 ± 0.012±g;g5o 
0.149 ± O.Olll^o 


3.4 • 10" 4 


7.3 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.140 ± 0.022lg;^| 
0.144 ±0.012l8;o°6 
0.136 ± O.OOSl^oos 


6.9 • 10" 4 


7.3 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.123 ± 0.02llg;^| 

0.122 ± o.oiiig;8^ 

0.141 ±0.008l8:8i4 


1.36 • 1(T 3 


7.3 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.158 ± 0.027l8;824 

0.140 ±0.01318:809 

0.149 ± 0.010l8:8o6 


2.67- 10" 3 


7.3 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.175 ± 0.070l8;?| 

0.193 ±o.o38l8:8n 

0.161 ± 0.02118:824 



Table 8: The leading proton production rate, r ' measured as a function of xl 
for protons with 0.04 < p\ < 0.15 GeV 2 , in bins of x and Q 2 , with averages (x) and 
(Q 2 )- Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Table 8 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


2.6 • 10" 4 


11 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.084 ± 0.021 1°;° 33 

o.isi ±0.019;^ 

0.140 ± 0.012±°;° 43 


4.6 • 10" 4 


11 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.134 ± 0.0201&009 
0.132 ± 0.01llg;oi6 
0.133 ± 0.007±jj;gjg 


9.2 • 10" 4 


11 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.129 ± 0.021±g;g?| 
0.140 ± O-OlllJooe 
0.150 ± 0.008±g;g^ 


1.83 • 10~ 3 


11 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.115 ± 0.018±g;ggf 
0.165 ±0.013±°;°°£ 
0.127 ± 0.007±g;gg| 


3.98 • 10~ 3 


11 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.110 ± 0.031±g;g?§ 
0.129 ± 0.018j£jj32 
0.157 ±0.0131^ 


5.1 • 10" 4 


22 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.112 ± 0.029ig:oo9 
0.133 ±0.0171^ 
0.148 ± 0.013±g;gg| 


9.2 • 10" 4 


22 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.146 ± 0.028±g;g|§ 

0.153 ±0.015ig:o25 
0.157 ±0.010±°;°° 4 


1.84 • 1(T 3 


22 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.101 ± o.oi9lg:8o6 

0.141 ± 0.013lo;o?2 
0.144 ±0.0091^ 


3.66 • 1(T 3 


22 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.177 ±0.036tJ];[! 4 8 
0.165 ±0.018^20 
0.156 ± O.Ollt^o 


7.83 • 10~ 3 


22 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.166 ± 0.046lg;oi4 
0.185 ± 0.031±g;g^ 
0.147 ±0.0151°;°°? 


1.03 • 10~ 3 


44 


0.38 - 0.56 
0.56 - 0.74 


0.128 ± 0.042lg;Soi 
0.090 ± 0.016lg;gg§ 
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Table 8 (cont.) 



(x) 


{Q ) {GeV z ) 


xl range 








0.74 - 0.92 


0.152 ±0.018j£|j42 


1.86 • 10~ 3 


44 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.190 ± 0.048±g;gf§ 
0.136 ±0.019j£roijj 
0.115 ± O-OlOloiool 


3.68 • 1(T 3 


44 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.082 ± 0.023±g;JJ§ 
0.156 ±0.023lg;^ 
0.158 ± 0.015lg;gS 


7.33 • 10~ 3 


44 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.138 ± 0.045±g;gJ5 

0.169 ±0.028±o;o22 
0.188 ± 0.020^056 


1.54- 1(T 2 


44 


0.38 - 0.56 
0.56 - 0.74 
0.74 - 0.92 


0.130 ± 0.055lg;g|g 
0.195 ±0.043^78 
0.118 ±0.016±g;gg 


2.00 • 10~ 3 


88 


0.38 - 0.65 
0.65 - 0.92 


0.092 ± 0.030±g;gf| 
0.109 ± 0.018^;^ 


3.59 • 10~ 3 


88 


0.38 - 0.65 
0.65 - 0.92 


0.236 ± 0.053±g;g|g 
0.153 ± 0.018±g;gJ| 


7.37 • 10" 3 


88 


0.38 - 0.65 
0.65 - 0.92 


0.190 ± 0.047±g;g|g 
0.191 ± 0.024±gio« 


1.42 • 10" 2 


88 


0.38 - 0.65 
0.65 - 0.92 


0.153 ± 0.040l{j;{]J! 
0.179 ± 0.026±g;gg| 


3.01 • 10~ 2 


88 


0.38 - 0.65 
0.65 - 0.92 


0.137 ±0.053lo;oio 
0.158 ± 0.032t°;° 3 ° 


4.U0 • 10 


237 


0.38 - 0.65 
0.65 - 0.92 


0.051 ± 0.048I qq 3 
0.067 ± 0.027±g:oM 


7.52 • 10~ 3 


237 


0.38 - 0.65 
0.65 - 0.92 


0.162 ± 0.053±g;g|| 
0.164 ± 0.027±g;g§e 


1.47- 10~ 2 


237 


0.38 - 0.65 
0.65 - 0.92 


0.161 ± 0.052±g;gg 
0.095 ± 0.015±g;gg| 


3.25 • 10~ 2 


237 


0.38 - 0.65 
0.65 - 0.92 


0.179 ±0.054i^ 
0.123 ± 0.019lg;g^ 
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(x) 


(Or) (GeV 2 ) 


xl range 




9.6 • 10~ 5 


4.2 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.114 ± 0.027±g:8i4 
0.116 db 0.033±g;g|| 
0.121 ±0.016±°;°f 4 


1.7- 10" 4 


4.2 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.102 ± 0.014±o;o26 
0.105 ± 0.0181&034 
0.133 ± 0.012t°;° 3 ? 


3.5 • 10" 4 


4.2 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.114 ±0.018+^2? 
0.212 ± 0.043l{j;?72 
0.099 ± O-OlOlgion 


6.9 • 10" 4 


4.2 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.097 ±0.01518$$ 
0.166 ±0.034ig;g|| 

0.151 ±0.016^20 


1.46 • 10~ 3 


4.2 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.106 ± 0.024+°;°^ 
0.090 ± 0.027lg;g|J 
0.154 ± 0.022lg;o33 


1.9 • 10" 4 


7.3 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.158 ± 0.034+°;o69 
0.068 ± 0.026±g;glg 
0.100 ± 0.017ig;^ 


3.4 • 10" 4 


7.3 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.138 ± 0.025±g;^| 
0.149 ±0.034lg:o?7 
0.110 ±0.013lg;oJ? 


6.9 • 10" 4 


7.3 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.137 ± 0.025±g;g^ 
0.098 ± 0.024±g;gg§ 
0.160 ± 0.0181^32 


1.36 • 1(T 3 


7.3 


0.62 - 0.72 
0.72 - 0.82 
0.82 - 0.92 


0.127 ± 0.024 jftooi 
0.170 ± 0.040±g;g|| 
0.127 ±0.016lg;g^ 


2.67- 10~ 3 


7.3 


0.62 - 0.77 
0.77 - 0.92 


0.245 ± 0.103±g:iao 
0.113 ±0.027lg;oi8 



Table 9: The leading proton production rate, r ^ >, measured as a function of xl 
for protons with 0.15 < p\ < 0.5 GeV 2 , in bins of x andQ 2 , with averages (x) and 
(Q 2 ) . Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Table 9 (cont.) 



(x) 


(Q ) (GeV 2 ) 


xl range 


r LP(3) 


2.6 


10" 


4 


11 


0.62 
0.72 


- 


0.72 
0.82 


0.170 ±0.047' 
0.088 ± 0.048" 


f0.063 
-0.036 
f0.225 
-0.012 










0.82 


: 


0.92 


0.169 ±0.027" 


f0.023 
-0.059 


4.6 


io- 


4 


11 


0.62 


- 


0.72 


0.176 ± 0.030 


f0.015 
-0.040 










0.72 


- 


0.82 


0.126 ±0.030" 


f0.076 
-0.074 










0.82 


- 


0.92 


0.134 ±0.015" 


f0.020 
-0.027 


9.2 


10" 


4 


11 


0.62 


- 


0.72 


0.142 ±0.023" 


f0.015 
-0.039 










0.72 


- 


0.82 


0.158 ±0.040" 


f0.026 
-0.117 










0.82 


- 


0.92 


0.148 ±0.016; 


fO.017 
-0.019 


1.83 


• 10 


-3 


11 


0.62 


- 


0.72 


0.107 ± 0.018±o:oi2 










0.72 
0.82 


- 
- 


0.82 
0.92 


0.071 ± 0.019" 
0.145 ±0.015" 


f0.089 
-0.008 
f0.020 
-0.017 


3.98 


• 10 


-3 


11 


0.62 


- 


0.77 


0.135 ±0.034 


f0.058 
-0.014 










0.77 


- 


0.92 


0.199 ± 0.033" 


f0.018 
-0.050 


5.1 


io- 


4 


22 


0.62 


- 


0.72 


0.125 ±0.033" 


fO.019 
-0.052 










0.72 


- 


0.82 


0.179 ±0.069; 


f0.038 
-0.121 










0.82 


- 


0.92 


0.127 ±0.020 


f0.053 
-0.012 


9.2 


io- 


4 


22 


0.62 


- 


0.72 


0.127 ± 0. 024" 


f0.057 
-0.016 










0.72 


- 


0.82 


0.106 ± 0.027" 


fO.014 
-0.019 










0.82 


- 


0.92 


0.147 ± 0. 018" 


fO.010 
-0.021 


1.84 


• 10 


-3 


22 


0.62 


- 


0.72 


0.084 ±0.017" 


f0.028 
-0.008 










0.72 


- 


0.82 


0.158 ±0.039" 


fO.017 
-0.063 










0.82 


- 


0.92 


0.157 ±0.021 


f0.016 
-0.036 


3.66 


• 10 


-3 


22 


0.62 




0.72 


0.115 ±0. 029" 


f0.019 
-0.052 










0.72 




0.82 


0.112 ±0.030" 


f0.097 
-0.011 










0.82 




0.92 


0.124 ±0.020" 


f0.018 
-0.068 


7.83 


• 10 


-3 


22 


0.62 




0.77 


0.110 ±0.037 


f0.064 
-0.024 










0.77 




0.92 


0.223 ± 0.044" 


f0.015 
-0.127 


1.03 


• 10 


-3 


44 


0.62 




0.77 


0.171 ± 0.054" 


f0.048 
-0.039 










0.77 




0.92 


0.100 ± 0.023" 


fO.014 
-0.024 


1.86 


• 10 


-3 


44 


0.62 




0.77 


0.124 ±0.030" 


f0.013 
-0.024 










0.77 




0.92 


0.141 ±0.022 


f0.068 
-0.015 
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Table 9 (cont.) 



(x) 


(Q 2 ) (GeV 2 ) 


xl range 


r LP(3) 


3.68 • 10" 3 


44 


0.62 - 0.77 
0.77 - 0.92 


0.170 ± 0.049lao58 
0.117 ±0.022±o;o?6 


7.33 • lO- 3 


44 


0.62 - 0.77 
0.77 - 0.92 


0.185 ± 0.078lg;J^ 
0.120 ± 0.024lg:o74 


1.54 • 10" 2 


44 


0.62 - 0.77 
0.77 - 0.92 


0.196 ± 0.078lg;rao 
0.182 ± 0.050lo;oro 


2.00 • 10" 3 


88 


0.62 - 0.77 
0.77 - 0.92 


0.156 ± 0.0771^8? 
0.125 ± 0.034±g;g^ 


3.59 • 10~ 3 


88 


0.62 - 0.77 
0.77 - 0.92 


0.208 ± 0.087lg;o44 
0.129 ± 0.027±g;gf| 


7.37 • 10~ 3 


88 


0.62 - 0.77 
0.77 - 0.92 


i n nriQ 

0.061 ± 0.031±g;g^ 
0.127 ±0.035lg;o4o 


1.4Z • 1U 


oo 


n £o c\ 11 

U.DZ - U. II 

0.77 - 0.92 


n i A k J- n nc;^+0.338 
U.140 ± U.UOD_q 

0.268 ± 0.079lg:264 


3.01 • 10~ 2 


88 


0.62 - 0.77 
0.77 - 0.92 


0.106 ± 0.063±g;g2o 
0.215 ± 0.1071^46 


4.00 • 10~ 3 


237 


0.62 - 0.92 


0.248 ±0.107lg;?^ 


7.52 • 10" 3 


237 


0.62 - 0.92 


0.183 ±0.049+^ 


1.47 • 10" 2 


237 


0.62 - 0.92 


0.203 ± 0.057lg;^ 


3.25 • 10" 2 


237 


0.62 - 0.92 


0.164 ±0.045i°;°55 
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(Q 2 ) (GeV 2 ) 


X 


r LP(2) 


4.2 


9.6 • 10" 5 
1.7- 10" 4 
3.5 • 10~ 4 
6.9 • 10" 4 
1.46 • 10~ 3 


0.256 ± 0.006lg;^ 
0.245 ± 0.004±g;g^f 
0.237 ± 0.004jlg;g$ 
0.234 ± 0.004+°;°^ 

0.214 ± 0.006±g;g^g 


7.3 


1.9 • 10" 4 
3.4 • 10" 4 
6.9 • 10" 4 
1.36 • 10~ 3 
2.67- 10~ 3 


0.246 ± 0.006±g:^2 
0.231 ± 0.005±g;gi§ 
0.236 ± 0.0051^2? 
0.243 ± 0.006±g;gi| 
0.252 ± 0.012+°;° 2 ^ 


11 


2.6 • 10" 4 
4.6 • 10" 4 
9.2 • 10" 4 
1.83 • 10~ 3 
3.98 • 1(T 3 


0.253 ± 0.007t|];j] 2 || 
0.241 ± 0.005lg;^ 
0.242 ± 0.005±g;gi| 
0.230 ± 0.005±g;gi| 
0.258 ± 0.008±g;g|g 


22 


5.1 • 10" 4 

9.2 • 10" 4 
1.84- 10" 4 
3.66 • 10~ 3 
7.83 ■ 10~ 3 


0.242 ± 0.007t[!;j] 2 3 
0.244 ± 0.006±g;gi| 
0.227 ± 0.005±g;g?| 
0.253 ± 0.007±g;g?g 
0.254 ± 0.010±g;g^ 


44 


1.03 • 10~ 3 
1.86 • 10~ 3 
3.68 • 10~ 3 
7.33 • 10" 3 
1.54- 10~ 3 


0.228 ± 0.009±g:gi| 
0.238 ± 0.008±g;gi| 
0.270 ± 0.009±g;g|? 
0.271 ± 0.011±g;gf? 
0.253 ± 0.014±°;°f^ 


88 


2.00 • 10~ 3 
3.59 • 10~ 3 
7.37 • 10" 3 
1.42 • 10~ 2 

3.01 • 10~ 2 


0.215 ± 0.013tg;gi9 
0.245 ± O.OlllJjjjg 
0.243 ± 0.013t[J;o 2 7 
0.265 ± 0.0151m 
0.241 ± 0.019±g;g|J 


237 


4.00 • 10" 2 
7.52 • 10~ 2 
1.47- 10~ 2 
3.25 • 10~ 2 


0.242 ± 0.0241^? 
0.264 ± 0.018±g;g§| 
0.227 ± 0.015±g;gf? 
0.232 ± 0.015+°;^ 



Table 10: The leading proton production rate, r LP ( 2 \ measured as a function of 
x in bins of Q 2 , with average (Q 2 ), for protons with 0.32 < xl < 0.92 and p\ < 0.5 
GeV 2 . Statistical uncertainties are listed first, followed by systematic uncertainties. 
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Q 2 (GeV 2 ) 


( r LP(2)} 




0.6 < x L < 0.97 


0.32 <x L < 0.92 


5.1 

15.8 

81.1 


0.145 ±0.0011^ 
0.149 ±0.0011°;^ 
0.151 ± 0.002+°;°^ 


0.238 ± 0.002+°;°^ 

0.241 ± 0.002jlg;gi5 
0.245 ± 0.003+°;°};* 



Table 11: The leading proton production rate, (r LF ^), averaged over x, as a 
function ofQ 2 , for protons withp\ < 0.5 GeV 2 in two Xl ranges as denoted in the 
table. Statistical uncertainties are listed first, followed by systematic uncertainties. 
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(Q 2 ) (GeV 2 ) 


X 


„LP(2) 
r 2 


4.2 


9.6 • 10" 5 
1.7- 10" 4 
3.5 • 10~ 4 
6.9 • 10" 4 
1.46 • 10" 3 


0.262 ± 0.006lg;^ 
0.234 ± 0.004l°;° 2 ^ 
0.195 db 0.003±g;gi| 
0.167 ± 0.003±g;g}§ 
0.133 ± 0.004±g;g^ 


7.3 


1.9 • 10" 4 
3.4 • 10" 4 
6.9 • 10" 4 
1.36 • 10~ 3 
2.67- 10~ 3 


0.290 ± 0.0071^ 

0.239 ± 0.005lg:o2o 
0.209 ± 0.004+°;°}^ 
0.183 ±0.004l^{ 4 
0.165 ± 0.008±g;gf| 


11 


2.6 • 10" 4 
4.6 • 10" 4 
9.2 • 10" 4 
1.83 • 10~ 3 
3.98 • 10~ 3 


0.327 ± 0.010±g;gP 
0.281 ± 0.005±g;gJ| 

0.239 ± 0.005lg;gi9 
0.191 ± 0.004+°;°}^ 
0.181 ± 0.006±g;gf| 


22 


5.1 • 10" 4 

9.2 • 10" 4 
1.84- 10~ 3 
3.66 • 10~ 3 
7.83 • 10" 3 


0.336 ± O.OlOl^f 3 , 
0.299 ± 0.007±g;g|| 
0.230 ± 0.005±g;g^| 
0.210 ± 0.006+^;^| 
0.176 ± 0.007^;^ 


44 


1.03 • 10" 4 
1.86 • 10~ 3 
3.68 • 10~ 3 
7.33 • 10~ 3 
1.54- 10~ 2 


0.322 ± 0.013±g:^f 
0.286 ± 0.009±g;g|§ 
0.265 ± 0.009j#Q22 
0.215 ± 0.0081^20 
0.164 ± 0.009±g;gy 


88 


2.00 • 10~ 3 
3.59 • 10~ 3 
7.37 • 10" 3 
1.42 • 10~ 2 

3.01 • 10~ 2 


0.286 ± 0.017lg;^? 
0.275 ± 0.013±g;g|g 
0.214 db 0.011±g;g|| 
0.190 ± 0.010±g;g|g 
0.138 ± 0.011±g;g^ 


273 


4.00 • 10~ 3 
7.52 • 10~ 3 
1.47- 10~ 2 
3.25 • 10~ 2 


0.287 ± 0.028lg;gp 
0.260 ± 0.0171^34 
0.173 ± 0.0121°;™ 
0.139 ± 0.0091°;™ 



Table 12: The leading proton structure function, F 2 , measured as a function 
of x in bins of Q 2 , with average (Q 2 ), for protons with 0.32 < xl < 0.92 and 
p\ < 0.5 GeV 2 . Statistical uncertainties are listed first, followed by systematic 
uncertainties. 
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Figure 2: Comparison between data (dots) and reweighted MC (shaded his- 
tograms) for DIS quantities for the LPS sample: (a) Z coordinate of the vertex, 
(b) E — Pz distribution, (c) energy of the scattered positron, E' e , (d) polar angle of 
the scattered positron, 6 e , (e) virtuality of the exchanged photon, Q 2 , (f) invariant 
mass of the hadronic system, W, and (g) Bjorken scaling variable, x. 
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ZEUS 




Figure 3: Comparison between data (dots) and reweighted MC (shaded his- 
tograms) for LPS specific quantities: (a) fractional longitudinal momentum, xl, 
(b) squared transverse momentum, p\, (c) minimum track distance from the edge 
of the pot, A pot , and (d) minimum track distance from the beampipe, A pipe . 
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Figure 4: Acceptance of the LPS spectrometers (a) sl23 and (b) s456 in the Xl, 
p% plane. The dashed lines delimit the three p\ integration ranges. 
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Figure 5: faj i?aiio pfjvc o/ number of events with a track in the LPS and a 
neutron candidate in the FNC to the number of events with a track in the LPS. The 
dots represent the data and the shaded histogram the reweighted MC; (b) expected 
fraction R of positive pions and kaons reconstructed in the LPS as a function of 
Xl- The error band reflects the statistical uncertainty derived from the LPS- FNC 
data sample. 



49 



ZEUS 




Figure 6: The double- differential cross-section cPctlp / ' dx^dp\ for Q 2 > 3 GeV 2 
and Ah < W < 225 GeV as a function of p\ in bins of Xl,. The circles and the dots 
are the ZEUS data measured with the spectrometers sl23 and s456, respectively. For 
clarity, only the statistical uncertainties are shown. The systematic uncertainties 
are listed in Table 1. The lines are the result of a fit to a function A ■ e~ h ' p T , 
as described in the text. The solid lines indicate the range in which the fit was 
performed. The bands show the statistical uncertainty of the fit. 
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Figure 7: The p\-slope, b, of the cross-section d? ! <7lp / ' dx hdp\ > as defined by 
the parameterisation A ■ e~ b ' PT and obtained from a fit to the data in bins of xl, 
in the kinematic range indicated in the figure. The bands represent the systematic 
uncertainty. 
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Figure 8: The p\-slope, b, of the cross-section cPctlp / 'dx^dp^, as defined by 
the parameterisation A ■ e~ b ' PT and obtained from a fit to the data in bins of xl, 
in different kinematic ranges as indicated in the figure. The bands represent the 
systematic uncertainty. For the ZEUS 1995 data [7] and the ZEUS 1994 data [16], 
the inner vertical bars represent the statistical uncertainties, the outer bars the 
statistical and systematic uncertainties added in quadrature. 
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Figure 9: The leading proton production rate, \/a- inc • aip/dxL, as a function of 
xl in three p\ ranges as indicated in the figure. The measurements in the higher 
p\ range are multiplied by a factor two for visibility. The bands represent the 
systematic uncertainty. 
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Figure 10: The leading proton production rate, \/a mc -dai,-p/dxL, fortwo ranges of 
Q 2 as indicated in the figure. The bands represent the systematic uncertainty. For 
the ZEUS 1995 data [7] the inner vertical bars represent the statistical uncertainties, 
the outer bars the statistical and systematic uncertainties added in quadrature. 
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Figure 11: The leading proton production rate, 1 / cr inc • da^p/dx^, for < 0.04 
GeV 2 in the kinematic ranges indicated in the figure. Other details as in Fig. 8. 
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Figure 12: The leading proton production rate, r LP ^ 3 \ as a function o/xl in bins 
of x and Q 2 , jor p\ < 0.5 GeV 2 . The bands represent the systematic uncertainty. 



56 



ZEUS 



1 1 

• ZEUS 12.8 pb" 1 






1 1 1 


1 


i i i 


i i i 
0.3 




i 








1 








f 


0.2 


0.32< x T <0.92 Q 2 =237 GeV 2 






















0.1 


p 2 < 0.5 GeV 2 












0.3 




• 






i 


+ : 
















0.2 


Q 2 =88 GeV 2 












0.1 












0.3 




♦ • 


» 




♦ 


r 


0.2 
0.1 




Q 2 =44 GeV 2 




















0.3 














• 


• 




♦ 


0.2 
0.1 






Q 2 =22 GeV 2 


















0.3 






















* 


• 




0.2 




















Q 2 =11 GeV 2 






0.1 












0.3 
























• 


♦ 


0.2 
0.1 










Q 2 =7.3 GeV 2 














0.3 












* 

• 


0.2 












Q 2 =4.2 GeV 2 ; 
















0.1 












1 1 






, , 1 









-4 -3 -2 

10 10 10 

x 

Figure 13: The leading proton production rate, r LP ( 2 \ as a function of x in 
bins Q 2 , in the kinematic range indicated in the figure. The bands represent the 
correlated systematic uncertainty. 
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Figure 14: The leading proton production rate averaged over x, (r LP ^), as a 
function of Q 2 in two different ranges of xl, in the kinematic range indicated in 
the figure. The bands represent the correlated systematic uncertainty. 
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Figure 15: The leading proton structure function, F LP ( 2 \ as a function of x 
in bins of Q 2 , in the kinematic range indicated in the figure. For clarity, only the 
statistical uncertainties are shown. The systematic uncertainties are listed in Ta- 
ble 12. The lines show the F 2 parameterisation scaled by 0.24, the approximate 
average value ofr LP ( 2 \ The bands show the one-standard deviation limits of the F 2 
parameterisation. 
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Figure 16: The rate l/cr inc ■ da^s/dxi for leading proton (dots) and leading 
neutron production (circles). The bands show the systematic uncertainties. 
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Figure 17: The p\-slope, b, as a function of xl for leading proton (dots) and 
leading neutron (circles) production. The bands show the systematic uncertainties. 
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Figure 18: A Regge-based model [13] compared to (a) the measured leading 
proton production rate, l/cr- mc ■ da^/dxi, and (b) the p^-slope, b. The bands show 
the systematic uncertainties. 
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Figure 19: Expectations of various Monte Carlo models of DIS, as described in 
the figure, compared to (a) the leading proton production rate, l/o"i nc ■ da^p/dxi, 
and (b) the p\-slope, b. The bands show the systematic uncertainties. 
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